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Abstract 

A system-reservoir nonlinear coupling model has been proposed for a situation where the reservoir is nonlinearly driven by an 
external Gaussian stationary noise which exposes the system particles to a nonequilibrium environment. Apart from the internal 
thermal noise, the thermodynamically open system encounters two other noises that are multiplicative in nature. Langevin equation 
derived from the resulting composite system contains the essential features of the interplay between these noise processes. Based on 
the numerical simulation of the full model potential, we show that one can recover the turnover features of the Kramers dynamics 
even when the reservoir is modulated nonlinearly by an external noise. 



1. Introduction 

Inspired by a work of Christiansen [Q]], where a chemical 
reaction was considered as a diffusion problem, Kramers J2l 
introduced a Brownian motion model in a one-dimensional 
(along the reaction coordinate) force field to predict the exis- 
tence of several kinetic regimes depending on the magnitude of 
the friction (very low or energy diffusion regime, and moder- 
ate to high or spatial diffusion regime). A clear understanding 
of the pre-factor of the Kramers equation is useful not only for 
completeness of the theory of escape rate, but also for explain- 
ing various phenomena. Therefore, for the last few decades, 
much effort has been put into extending the Kramers model 
Many authors have devised methods for 
obtaining escape rate in the whole range of friction by ex- 
tending the basic assumptions found in the original Kramers 
work, known as the Kramers turnover problem |3|5l Tol IT | . 
Kramers has shown that the rate constant is proportional to y 
(dissipation constant) when y is low and proportional to y~ l 
when y is high. It can thus be expected that the value of 
the rate constant reaches a maximum at an intermediate y and 
decreases to zero when y approaches either zero or infinity. 
This dependence of the rate constant on friction is known as 
Kramers turnover. However, Kramers could not derive a uni- 
form expression for the rate, valid for all values of the friction 
coefficient. In fact, analytical solutions for Kramers equation 
are only possible for very simple interaction models Q3], |4j |6[]. 
The systematic solution of the Kramers turnover problem for 
the thermodynamically closed system was given by Pollak- 
Grabert-Hanggi(PGH) lfl2Tl (one of the foremost studies about 



turnover) who generalized the Kramers model to an arbitrary 
time-dependent friction and demonstrated that the turnover for- 
mula due to Mel'nikov and Meshkov can be obtained 
without any ad hoc bridging. Later, the PGH theory was gener- 
alized to many dimensions jl4[l5ll . 

Activated rate processes in one-dimensional surface diffusion 
have been studied by Pollak et al. Il4 llol . Hershkovitz and 
Pollak 11 51 studied the length dependence of the classical acti- 
vated transfer rate across a bridge and found that the Kramers 
turnover theory in the rate suffices for understanding the bridge 
length and friction dependence of the rate. Rips and Pollak 11711 
extended the PGH method in the context of quantum Kramers 
turnover problem. Segal et al. 11811 have provided the first anal- 
ysis of the transition from the tunnelling to the thermally acti- 
vated regime in a variant of the quantum Kramers problem as 
a function of the barrier length. Vega et al. lfl9ll have stud- 
ied the Kramers turnover theory in activated atom-surface dif- 
fusion using mean first passage time. Shepherd and Hernan- 
dez I 20l 21. 221 have exploited the mean first passage time 
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(MFPT) based rate formula to analyze the interplay between 
Kramers turnover and resonant activation for the escape rates 
on stochastic bistable and aperiodic potentials. The develop- 
ment in Ref. lE-ll is particularly interesting as it investigates the 
low friction regime (the most difficult part of Kramers turnover 
theory to illustrate) in conjunction with stochastic aperiodic po- 
tentials. Recently, Kramers turnover has also been realized dur- 
ing the investigation of the forward and backward reaction rates 
of the LiNC^LiCN isomerization reaction in a bath of argon 
atoms at various densities using molecular dynamics simula- 
tions due to Garcfa-Muller et al. 12311 . Their work provides 
clear evidence for the increase in rates with microscopic fric- 
tion in the energy-diffusion regime in chemical system. 

The last few decades have observed a crescendo of research 
activity in the field of nonequilibrium statistical mechanics us- 
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ing the system-reservoir (SR) model 11241 |25[. In the over- 
whelming majority of situations, the interaction between the 
system and the reservoir has been considered to be linear in 
bath co-ordinates as well system co-ordinates. This in turn 
relates the additive noise of the thermal bath with linear dis- 
sipation of the system through fluctuation-dissipation relation 
(FDR). On the other hand, if the SR coupling is nonlinear in 
system coordinate, the corresponding Hamiltonian gives rise to 
a Langevin equation with state-dependent dissipation and in- 
ternal multiplicative thermal noise. As the total SR combina- 
tion is thermodynamically closed, the energy balance condition 
again is reflected through FDR II26I1 . If density of bath modes 
is such that the associated noise is stationary and Gaussian, one 
can numerically solve the associated Langevin equation for bar- 
rier crossing dynamics to observe the turnover phenomena. It 
should be recognized that it is very hard to obtain a simple ex- 
pression for escape rate, even for a white noise process, when 
the dissipation is state dependent and the noise process is mul- 
tiplicative in nature. Consequently, PGH II 1211 -type analysis for 
turnover problem is very hard to achieve in such cases 023J. At 
this point, one might wonder if the turnover phenomena can be 
observed when the SR combination is thermodynamically open 
and hence there is no energy-balance relation like FDR. 

Among many other situations, the SR combination will be 
thermodynamically open if one drives the system externally 
(keeping the reservoir in thermal equilibrium). On the other 
hand, in spite of direct driving, one may expose the reservoir 
to an external modulation. A number of different situations de- 
picting the modulation of the bath by an external noise may be 
physically relevant |27, 28, 29[ 3(i • Whether the system or the 
reservoir is driven by an external noise, there is an additional 
mechanism to inject energy into the system and clearly, there 
is no FDR in such a situation. Inspection of any such situa- 
tion may be relevant to examine the turnover phenomena in the 
rate. In what follows, we address the later situation where the 
reservoir is modulated externally by a random force to make 
the SR combination thermodynamically open. It is our aim 
here to search for the signature of Kramers turnover in the rate, 
emerging from the nonlinear driving of the bath by an exter- 
nal noise. The effect of nonlinear modulation of the reservoir 
by an external agency is considered indispensable in explaining 
the phenomena of activation of a quasibound species (reactants 
surrounded by the solvent molecules) above its trapping poten- 
tial in the presence of high intensity light sources. For example, 
one may consider an isomerization reaction (A ^ B) in a pho- 
tochemically active solvent in the presence of an external light 
source with high intensity. 

2. Methodology 

To start with, we consider a classical particle of unit mass 
being coupled to a heat bath consisting of N-mass weighted 
harmonic oscillators, characterized by the frequency set {a>j} 
(i.e. the bath degrees of freedom are described by an ensemble 
of oscillators). In addition to that, the heat bath is nonlinearly 
driven by an external noise identified as e(t). The Hamiltonian 



for the composite system is 

H — Hs + Hb + Hsb + Hi nt 



(1) 



where the Hamiltonian of the system is expressed as: Hs = 
(p 2 /2) + V(q) with V(q) being the potential energy function and 
p and q being respectively the coordinate and the momentum of 

the system particle. H b +Hsb = Yj% 



) [xj-Cjf(q)} 



where {xj,pj} are the variables for the j-th bath oscillator. 
The system-heat bath interaction is given by the coupling 
term CjU)jf{q) where Cj is the coupling strength and f(q) is 
some well-behaved function of the system coordinate q only. 
Through the insertion of the term f(q), we have considered the 
SR interaction to be, in general, nonlinear. For bilinear system- 
bath coupling, f(q) would have been taken as some linear func- 
tion of q. The interaction between the heat bath and external 
noise e(t) is taken as Hi nt = K j8 ( x j) e (0 where Kj denotes 
the strength of the interaction and gixj) is an arbitrary analytic 
function of the bath variable Xj. This type of interaction makes 
the bath variables explicitly time dependent. A large class of 
phenomenologically modelled stochastic differential equations 
may be obtained from a microscopic Hamiltonian for a partic- 
ular choice of coupling function g{xj) and have already been 
used for microscopic realization of Kubo-type oscillator and 
correlated noise processes 32 1. If one chooses, for exam- 
ple, g(xj) = jxj, the spring constants of the bath oscillators be- 
come fluctuating. In what follows, we choose g(xj) = xj + jx 2 ; 
a linear-linear(LL) and a square-linear(SL) coupling between 
the noise and bath variables. Recently, the effect of such LL 
coupling and SL coupling between the system and the reser- 
voir have been studied by Tanimura and coworkers rf33l in the 
context of spectroscopic studies. The notation of the hierar- 
chical representation of the underlying Hamiltonian system in 
Eq.(Q~|i and its subsequent analytic representations in extended 
Langevin equations has been discussed by Popov and Hernan- 
dez H- 

In what follows the external noise e(t) is taken to be station- 
ary, Gaussian with statistical properties (e(t)) e = 

{e(t)e(t')) e = «A(f - = ^ exp [~-^ . ( 2 ) 

where D £ is the strength of the noise and T e is its correlation 
time. For t £ — > 0, e(f) becomes <5-correlated with statistical 
property: {e(t)e(t')) e - 2D e 5(t - t'). In Eq.©, <■ ■ ■ ) e implies 
averaging over the external noise processes. In this context, we 
want to mention that the presence of noise and nonlinearity are 
unavoidable in general physical systems, So, one must take into 
account the interplay between these two factors on the dynam- 
ics of the system. 

From Eq. ([TJ, we have the dynamical equations for the sys- 
tem and the bath variables as 

q(t) = -V'(q(t)) + f'(q(t)) £ Cj Cj) (x/t) - c jf(q(t))} , (3) 

j 

Xj(t) + [co 2 j + Kjeit)} Xj (t) = -Kje(t) + CjL0 2 jf(q(t)). (4) 
To solve Eq.© for Xj(t), we assume a solution of the form 

Xj(t) = X )(t) + Kj X )(t), (5) 



where x°(f) is the solution of the unperturbed equation of mo- 
tion (EOM) 



x)(t) + oj)x)(t) = cjcop(q(t)). 



(6) 



We now consider that at t — 0, the heat bath is in thermal equi- 
librium in the presence of the system but in the absence of the 
external noise e{t). Subsequently, at t — + , the external noise 
agency is switched on and the heat bath is modulated by e(t). 
Then x\(t) must satisfy the equation 



x)(t) + cj)x)(t) = -6(0 - e(t)x){t). 



(7) 



with the initial conditions xj(0) = p](0) = 0. Now, the solution 
of Eq.Q is given by 



x){t) 



— f dt' sinaj ,(t - t')e(t') 
Uj Jo 

1 f 



i r 

dt sin aidt- t')x%t)e(t'). 

"J Jo 1 

The formal solution of Eq.© is given by 

p°(0) 

x°.(t) = x l '(0) cos oj Kf) + — — sin to dt) 

J J J J 



(8) 



+CjO)j j dt' sin a>j 
Jo 



(t-t')f(q(t')), 



(9) 



where x°.(Q) and p j(Q) are respectively the initial position and 
momentum of the j-th bath oscillator. Now, using this solution 
in Eq.©, we have (after an integration by parts) the EOM for 
bath variables xj(f) [from Eq.©] as 

Xj(t) - Cjf(q(t)) = {x°j(0) - cjf(q(0))} cos u>j(t) 
P°j(0) 



-cj f dt' cos a)j(t - t') f'(q(t'))q(f) 
Jo 

k r' 

— - dt' sin a>i(t- f')e(f') 
Jo 

-— I dt' sin cjjit - f')^ l '(f')e(f' 
tt] Jo ; 



). 

(10) 

Using the above solution in Eq.Q, we obtain the EOM for the 
system variables as 

q(t) = -V'(q(t))-f'(q(t)) f dt'y(t - t')f'(q(t'))q(t') 

Jo 

+f'(q(t))F(t) + f'{q(t))n{t) - f(q(t)) 



: J dt' |^ cjKjOjj sin cjj{t - f')^(f')| e(f), 



(ID 

where the damping kernel is given by y{t) = £ j Cj(o 2 . cos a>jt. 
F(t) is the internal thermal noise generated through the coupling 



between the system and the heat bath and is given by 

F(t) = J^CjCojUx^m-CjfiqiO^cosajjt 

j 

P°j(0) 



(12) 



and 



n(t) = f dt'ip(t - t')e(t') 
Jo 

dressed noise that depend: 
ip(t) - ^ c jKjWj sin a>jt. 



')■ (13) 
is a dressed noise that depends on the external noise e(f) and 



(14) 



Clearly, the system does not encounter the external noise e(t) 
directly, rather, the driving of the bath by the external noise e(f) 
results in a dressed noise. The form of Eq. dTTl therefore sug- 
gests that the system is driven by two forcing functions F(f) and 
n{t). F(f) depends on the initial conditions of the bath oscillators 
for a fixed choice of the initial condition of the system degrees 
of freedom. To define the statistical properties of F(f), we as- 
sume that the initial distribution is the one in which the bath is 
equilibrated at t = in the presence of the system but in the 
absence of the external noise agency. Let us now digress a little 
bit about n(t). The statistical properties of n(t) are determined 
by the normal-mode density of the bath frequencies, the cou- 
pling of the system with the bath, the coupling of the bath with 
the external noise, and the external noise itself. Equation ( fT3l 
is reminiscent of the familiar linear relation between the polar- 
ization and the external field, where n and e play the role of the 
former and the latter, respectively. The function <p(f), thus may 
be taken as the response function of the bath. The very structure 
of n(t) suggests that this forcing function, although originating 
from an external force, is different from a direct driving force 
acting on the system. The distinction lies at the very nature of 
the bath characteristics (rather than system characteristics) as 
reflected in the relations Eqs. (fTJt and (fT4l . At this point, we 
note that the forcing term F (t) is deterministic. It ceases to be 
deterministic if it is not possible to specify all the jk;"(0)'s and 
p'!(0)'s, i.e., the initial conditions of all the bath variables, ex- 
actly. The standard procedure to overcome this difficulty is to 
consider a distribution of x°(0) and p°j(0) to specify the statis- 
tical properties of the bath-dependent forcing term F(t). The 
distribution of the bath oscillators is assumed to be a canonical 
distribution of the Gaussian form 



Peq (0) 



A^exp 
1 



fpf(0) 



+ -OJ 2 j (x° j (0)-C j f(q(0))y 



(15) 



where is the normalization constant. This choice of the dis- 
tribution function of the bath variables makes the initial noise 
Gaussian. It is now easy to verify the statistical properties of 
F(t) as (F(t)} = and (F(t)F(t')} = k B Ty(t - ?') where k B 



is the Boltzmann constant and T is the equilibrium tempera- 
ture. Here, (• • • } implies the average over the initial distribution 
given in Eq.(fT3b. The second relation is the FDR II 1 CHI which 
ensures that the bath was in thermal equilibrium at t = 0, in 
presence of the system. To proceed further, we consider the last 
term of Eq. dTTb as 

T(t) = f'(q(t)) f dt' V CjKjUj sin a)j(t - t')e(t')x Q (t'). (16) 
Jo } 

We now put the expression for x°j(t') from Eq.®. The solu- 
tion Eq.©, consists of two parts, the homogeneous solution of 
Eq.© which is the free evolution of bath variables is the fast 
part. The second one is the solution of the corresponding inho- 
mogeneous equation which gives the the forced oscillation ex- 
pressed as CjOjj L dt' sin a>j(t - t')f(q(t')). As the fast part dies 
out quickly for damped driven oscillator, we pick the particu- 
lar solution of Eq.© only for jc°(f) and consequently, Eq. dTTI) 
becomes 

q(t) = -V\q(t))-f(q(t)) f dt'y(t - t')f'(q(t'))q(t') 

Jo 

+f'(q(t))F{t) + f'(q(t))n(t) 

-f'(q(t)) f dt'e(t') f dt"f{q{t")) 
Jo Jo 

X |^ c2 j K j aj2 j sm w / f _ O sm w ;'( f ' _ -(17) 

To identify Eq.dTTb as a generalized Langevin equation, we 
must impose some conditions on the coupling coefficients Cj 
and Kj, on the bath frequencies uj and on the number N of the 
bath oscillators that will ensure that y(t) is indeed dissipative 
and the last term in Eq.(fT7b is finite for N — > oo. A sufficient 
condition for y(f) to be dissipative is that it is positive-definite 
and decreases monotonically with time. These conditions are 
achieved if N — > oo and if CjOpj and coj are sufficiently smooth 
functions of / l35Tl . As N — > oo, one replaces the sum by an 
integral over to weighted by a density of states p{to). Thus, to 
obtain a finite result in the continuum limit, the coupling func- 
tion Cj = c(a>) and kj = k(co) are chosen as c(a>) = and Similarly, 
k((jS) = ko where cq and Kq are constants and t c is the correlation 
time of the heat bath. The choice k{o>) = Kq is the simplest one 
where we assume that every bath mode is excited with the same 
intensity. This simple choice makes the relevant term finite for 
jV — > oo. Consequently, y(t) becomes 



7(f) = 



( r 2 \ 



da>p(tL>) cos u>t, 



(18) 



where l/r e may be characterized as the cut-off frequency of 
bath oscillators. The density of modes of p{to) of the heat bath 
is assumed to be Lorentzian, 



(19) 



This type of choice of p(w) may be encountered in many sit- 
uations in chemical physics and condensed matter physics[36, 



(;) 


r o? i 




T~ 2 + CO 2 



37L b&l I39TI and resembles broadly, in behavior, the hydrody- 
namic modes in certain macroscopic svstems lHOIl . With these 
forms of p(w), c(oj) and k(cj), we have the expression for y(t) 

as y{t) = y exp (—7-) which reduces to y{t) = 2c^5{t) = 2y5{t) 
for r e — » where y — and is a Markovian dissipation con- 
stant and consequently, one obtains <5-correlated internal noise 
processes. With these forms of density of modes p(a>) and cou- 
pling functions, c(co) and k((l>), the response function (p(t) can 
be written in the continuum limit as 

<p(t) = I da>p( <x>)c((jj)k{<x>)(jl> sin a>t 
Jo 



2 1 r°° sine 
— CqKq — I da>a> — - — 
t t c Jo T7 + 



a>t 

u 2 



C K I -t 

■exp — 



(20) 



Clearly, for t c — > 0, <p(f) reduces to ip(t) — 2cQKod(t). Now, us- 
ing the standard trigonometric identity, the last term in Eq.(IT7t 
can be written as 



A(0 = f'(q(t)) 
X 



i f dt'e(t') f dt"f[q(t' 
z Jo Jo 

2^ c 2 jKjD 2 j cos a>j(t - 2f + t") 

j 

-i f dt'e(t') f dt"f[q(t")] 
1 Jo Jo 

x ^ c ' 1 j K i CL>1 j cos w j( f _ O 



)] 



(21) 



Now, using the assumed expressions for the coupling functions 
c(a>) and k(u>) and the density of modes picS), one easily ob- 
serves that the two sums in Eq.(f2TJl may be approximated as a 
(5-function, 



} t c 2 jKjtu 2 j cos LOj(t + t" - 2t') = I dtop(u)) {c(io)} 2 k(lo) 

j ' J 

Xw 2 cos a)(t + t" - 2t') 



2c 2 K 6(t + t" - 2f')(22) 



^ c-jKjOj 2 cos ojjit - t") = 2clK 6(t - t"). (23) 
Thus, in the continuum limit, the expression for A(f) reduces to 

A(f) = c 2 K f'(q(t)) \ f dt'e(t') f dt"f[q(t")} 
Jo Jo 



x6(t + t" - 2t') 

- f dt'e(t') f dt" f[q(t")]5(t - t") 
Jo Jo 



(24) 



With the property of (5-function, the first double integral in 
Eq.(t24b may be written as 

f dt'e(t') f dt"f[q(t")]5(t+t"-2t') 
Jo Jo 



f(q(y))- 



(25) 





Figure 1: Domain of integration of Ldlfe(f) L dt" f[q{t")]S{t - f) in Eq. 

E). 

As the system variable evolves much slowly in comparison to 
the external noise e(f), the right hand side of Eq.(l25l> may be 
approximated as ^ |^'flfye(^)]/(g(0)). For large /, we note 
that as 

lim - [ dt'e(t') = (e(f)) e = 0, (26) 
r-xx, t Jo 

the first term in the expression of A(f) vanishes. 

To perform the second integration; 

j^dt'e(t')j^ dt" f[q(t")]6(t - t'), we consider the region 
of integration, shown as the shaded triangle in Figure 1. 
From the property of (^-function, one observes that the above 
integral will contribute only when t" = t but the inequality 
< t" < t' < t demands that at the same time, t' should be 
equal to t. Thus, the contribution from the integral come out 
only at point P and the value of this contribution is f(q{t))e(t). 
Using all the above facts, we obtain from Eq.dTTb the EOM for 
system variables, in the limit r e — > 0, as 

q(t) = -V'(q(t))-y[f'(q)] 2 p+f'(q(t))F(t) 

+f'(q(t))n(t) + JK f(q)f'(q)e(t). (27) 

This equation can be used to explore the distinctive aspects of 
the reservoir (driven nonlinearly by an external noise) modu- 
lated dynamics of the system in contrast to direct driving of the 
system by the external noise. This will help us to elucidate the 
special role of the reservoir response function in controlling the 
escape of a Brownian particle from the metastable state. 

3. Results and Discussion: Kramers turnover 

Before examining the noise induced transport, it is instruc- 
tive here to have a close look at the above Langevin equation, 
where three noise processes appear and all these noise pro- 
cesses appear multiplicatively. F(t) is the internal thermal noise 
for which FDR exists. n(t) is the dressed noise and e(f) is the 
external noise. Instead of nonlinear SR coupling, if one consid- 
ers bilinear coupling, i.e., f(q) - q, the above equation Eq.(l27l> 
reads as 

q(t) = -V'{q{ij) - yp + F(t) + Ml) + yme{t), (28) 



which indicates that both the thermal noise and dressed noise 
appear additively but the last noise containing term appears 
multiplicatively. The effect of interference of colored additive 
and multiplicative white noises on escape rate has also been ex- 
plored using this type of equation[4l]]. Let us now discuss a lit- 
tle bit on the origin of the noises appeared in Eq.(l28V F(t), the 
usual thermal noise appears due to the system-bath interaction. 
The driving of the reservoir by external noise yields the last two 
terms in Eq.(l28li. If we choose the bath-noise coupling func- 
tion g(xj) to be linear in bath variable, one will encounter the 
n(t) noise in Eq.(l28ll only and the last term will disappear. On 
the other hand, if g(xj) be quadratic, i.e., g(xj) = (l/2)x 2 j, n(i) 
term disappears and the last term plays its role in the dynamics. 
Here, it is interesting to note that the multiplicative nature of 
the last noise process stems from the nonlinear driving of the 
bath but not from the nonlinearity of system-reservoir coupling 
function, which is the case for the other two noises. Here, we 
enunciate a system without proof that if the bath-noise coupling 
function be g(xj) - axj + bx 1 . + cx^. + then for linear system 
reservoir coupling, the resulting Langevin equation will read as 

q(t) = -V'(q(t)) - yp + F(t) + n(t) + Bqe(t) + Cq 2 e(t). (29) 

At this point, it is instructive to consider the statistical prop- 
erty of the dressed noise n(t) which can be easily verified as 
(n(t)) e = and 

(n(f)7r(t')) e = f dt" f dt"'<p(t-t")<p(t' -f") 
Jo Jo 

xiK?" - *"')• (30) 

If we assume that the external noise e(f) is <5-correlated, i.e., 
<e(f)e(f')> = 2D e 6(t - t'), then in the limit t c -> 0, the corre- 
lation function of n(t) becomes (w(f)w(f')) c = 2yK^D e 6(t — t'). 
In passing, we observe that the system encounters an effective 
Gaussian additive noise £(f) [= F(t) + n(tj\ and another noise 
which appears multiplicatively. The noises %(t) and e(t) are not 
statistically independent, their correlation may be expressed as 
<£(f)e(f')> = <e(0£(O> = P(t - f) which one may calculate for a 
particular ifr(t). Thus, the two mutually correlated noises appear 
in the dynamical equation of the open system. The appearance 
of cross-correlated noises has already been encountered while 
explaining various physical phenomena BOl l42l |43j |44|]. Now, 
in terms of an auxiliary function G(q) and a Gaussian stationary 
noise R(t), the Langevin equation Eq.(l28l can be written as 

q(t) = -Y(q(t)) - yp + G(q)R(t), (31) 

with 

((R(f)) e = ((R(t)R(t')}} = 2S(t - t'), (32) 

where ((• • ■ )} implies average over the noise process /?(f)(this 
averaging over R(t) consists of two independent averaging, one 
over thermal noise F(T) and another over external noise e{t)). 
In Eq.lEU, 

G(q) = [(yk B T + y«%D e ) + D^^q 1 + 2yK () D e q\' 1 . (33) 

Clearly, Eq.(l30l> along with Eq.OTb. is not the FDR but serves 
as the thermodynamic consistency relation. 
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Figure 2: Change of the transition rate, k (sec ') with the dissipation constant, 
y (sec -1 ) for various temperatures in conjunction with ks=l, A=3.0, B =0.1 
k =0.05, and D e =5.0. 



recovery of Kramers turnover for an thermodynamically open 
system is the key issue of our present investigation. Figure 2 
also shows that for a given value of y, the escape rate increases 
with increase in the temperature, as it should be. With increas- 
ing temperature, the sharpness of the turnover of the escape rate 
also increases. 

To this end we would like to mention the works of Zhou [461 



and Kalmykov et al.\Al\. In both of the works, the authors have 
considered the standard Langevin equation with constant and 
additive 5-correlated white noise which relates with the dissipa- 
tion by means of FDR (and hence describe thermodynamically 
closed system). In the work of Zhou||46J, the Langevin equa- 
tion was solved numerically to study the nature of the barrier 
dynamics, whereas the matrix-continued fraction method has 
been exploited to examine the thermally activated escape from a 
double-well potential for all values of dissipation by Kalmykov 
et al. 14711 . In both the works, inevitable Kramers turnover was 
examined and compared with those obtained by the Mel'nikov 
and Meshkov method ll48ll . On the other hand, our present work 
deals with Kramers turnover in the case of open system in con- 
junction with both additive and multiplicative noises. 



We now proceed to examine the noise induced transport. To 
do this, we numerically solve the Langevin equation, Eq.(l28l) 
[considering only quadratic bath-noise coupling, g(xj) = bx 2 j], 
by the method developed by Sancho et al. for multiplica- 
tive noise and routinely calculate the MFPT ll45ll . the inverse of 
which gives the escape rate from the metastable potential well. 
In our numerical implementation, we consider a double well 
potential of the form: 

A , B a 
V(q) = --q +4?' -oo<4<+oo 

In Figure 2, we have plotted the rate k, obtained from 
Langevin simulation using the concept of mean first passage 
time, as a function of dissipation constant(y) for various tem- 
peratures. For small y, we observe that the rate increases with 
increase in y whereas, for moderate to large y, k decreases: 
the rates turnover with the (microscopic) friction (Figure 2). 
This observation can be explained with the help of the fact that 
the interaction between the system (say reactants) and the bath 
must transfer sufficient energy to activate the reactants above 
the energy barrier leading to products. The corresponding rate 
should therefore increase with the coupling represented by fric- 
tion. An increase in friction, however, also slows down the re- 
actants and induces a competing mechanism that reduces the 
rate. Thus the topology of the variation of k with dissipa- 
tion constant in the present work also exhibits a typical sig- 
nature of Kramers turnover. It is thus important to note that 
the simulation of the barrier crossing dynamics of the external 
noise-driven-reservoir-modulated dynamics of the system cap- 
tures the essential turnover features of the Kramers dynamics 
of the closed system. In the detailed balance principle, when 
instead of additive internal thermal noise (for which FDR ex- 
ists), the system encounters another multiplicative nonthermal 
noise that originates due to the modulation of the bath by an ex- 
ternal noise, one recovers Kramers turnover nature. Thus, the 



4. Summarizing Remarks 

Many physical processes (with arbitrary complexity) influ- 
enced by the surroundings can be modeled as a potential barrier 
crossing event. Kramers showed that there is a qualitative dif- 
ference in the barrier crossing dynamics at the low and high 
friction limits. Many authors have devised theoretical and com- 
putational models to describe the Kramers turnover by extend- 
ing the basic assumptions found in the original Kramers work. 
The open question to be addressed here is whether the Kramers 
turnover is realizable in that class of thermodynamically open 
systems when the reservoir is modulated nonlinearly by an ex- 
ternal noise and hence is relevant to chemical dynamics, in con- 
junction with other physical processes. 

This work is a continuation of our studies on the models to 
describe the Kramers turnover. In Ref. [49], Ray Chaudhuri et 
al. shown numerically that the well known Kramers turnover 
phenomena is restored when the bath is linearly modulated by 
an external noise. However, in the present work the bath is be- 
ing driven nonlinearly by an external noise. In this case, in spite 
of having a linear system-bath interaction, the nonthermal noise 
will appear multiplicatively in the Langevin equation. The ori- 
gin of this multiplicative nature lies in the nonlinear driving of 
the bath itself. We have also envisaged the Kramers turnover 
phenomenon for the present model. Main results of this work 
are presented in Figure 2 which show the behavior of the rate 
constant as functions of the friction coefficient of the environ- 
ment. From the aforesaid, we are led to the conclusion that ir- 
respective of the mode in which the bath nonequilibration takes 
place, the turnover phenomenon will make its appearance, and 
it is not only the additive noise that leads to such an observation, 
but also the multiplicative noise too has the potential to induce 
turnover. The observations of the present work are valid for all 
types of processes in which a classical system in contact with 
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a thermal heat bath is driven out of equilibrium by classical, 
generally time-dependent fluctuating forces. 
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